Tensor Rank and Stochastic Entanglement Catalysis for Multipartite Pure States 
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The tensor rank (also known as generalized Schmidt rank) of multipartite pure states plays an 
important role in the study of entanglement classifications and transformations. We employ powerful 
tools from the theory of homogeneous polynomials to investigate the tensor rank of symmetric states 
such as the tripartite state = --i=(|100) + |010) + |001)) and its A^'-partite generahzation |V[^]v). 
Previous tensor rank estimates are dramatically improved and we show that (i) three copies of IW3) 
has rank either 15 or 16, (ii) two copies of |Wjv) has rank 3A^ — 2, and (iii) n copies of \Wn) has 
rank 0{N). A remarkable consequence of these results is that certain multipartite transformations, 
impossible even probabilistically, can become possible when performed in multiple copy bunches or 
when assisted by some catalyzing state. This effect is impossible for bipartite pure states. 

PACS numbers: 03.67.-a, 03.65. Ud, 03.67.Hk 



Multipartite entanglement has attracted increasing at- 
tention due to its intriguing properties and potential ap- 
plications in both quantum information processing and 
condensed matter physics [H, . A central question in 
the subject concerns the convertibility between different 
multipartite entangled states by using local operations 
and classical communications (LOCC). If such a protocol 
is only stochastic (i.e., occurs with a non-zero probabil- 
ity) then we say that the two states are convertible via 
stochastic LOCC (SLOCC); when the transformation is 
reversible, the two states are called SLOCC equivalent. In 
bipartite systems, SLOCC convertibility is characterized 
by the Schmidt rank of the state: bipartite |?/') is SLOCC 
convertible to |(/)) if and only if the Schmidt rank of 
is no smaller than that of {(p). 

A generalization of the Schmidt rank in multipartite 
systems and also relevant to SLOCC transformations 
is the tensor rank. Formally, for states in A^-partite 
quantum systems, each of which is described by a d- 
dimcnsional Hilbcrt space Hi {i — 1,...,N), the ten- 
sor rank i-k{t(j) of a state Ii/j) G {S>q=i ^q, defined as the 
smallest number of product states {0^^^ \(t>f)}i=i...rk(ii) 
whose linear span contains \ip). The tensor rank has been 
extensively studied in algebraic complexity theory , 
and while it is easy to compute for = 2 (Schmidt rank), 
even for A'' = 3, determining the rank of a state is NP- 
hard [H. This is one reason why SLOCC convertibility 
in multipartite systems is so challenging. 

Despite this general difficulty, the analysis becomes less 
formidable when certain classes of states are considered 
such as symmetric states, i.e. those invariant under any 
permutation of its parties. Recently, symmetric states 



have received much attention in the study of entangle- 
ment measures 0, 0] and bound entanglement Q . Fur- 
thermore, the entanglement transformationproperties of 
symmetric states have been investigated |9| and exper- 
imental procedures have been designed which use sym- 
metric states in generating families of multiqubit SLOCC 
equivalent states [3 El- 



There is a natural correspondence between symmet- 
ric tensors and symmetric polynomials where the the- 
ory of homogeneous polynomials can be used to study 
the latter. As we will explore in greater detail, every 
homogenous polynomial possesses a quantity called the 
polynomial rank which is closely related to the tensor 
rank. The relationship between the two ranks allows for 
known results on the polynomial rank to be used directly 
on tensor rank estimations [l2l |. This method will prove 
to be quite powerful. 

It is easy to see that the tensor rank is an SLOCC 
monotone: if It/:) can be transformed into via SLOCC, 
then rk('0) > rk ((/>). In general the converse is not 
true however for any state SLOCC equivalent to the 
d-level iV-partitc GHZ state \GHZ%) = ^ E^i K)^^. 
tensor rank does decide convertibility [l^ . 

Observation 1 A CHZ-equivalent state \ipGHz) can be 
SLOCC transformed into iff vk^TpGHz) > i'k((/>). □ 

Two related types of phenomena studied in entangle- 
ment theory are multi-copy and entanglement-assisted 
entanglement transformations. Given a source state 
I?/;) and a target state {(p), if there is an integer k 
such that the transformation of Ith)'^'' to IS)^'' can be 
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achieved by LOCC, then we say that can be trans- 
formed to by multiple- copy entanglement transforma- 
tion (MLOCC). Similarly, if there is a state |c) such that 
the transformation of |'0) (S) |c) to \(j>) ® \c) is possible 
by LOCC, then we say that can be transformed to 
|(/)) by entanglement- assisted (or catalytic) transforma- 
tion (ELOCC). The state |c) is called a catalyst for the 
transformation. For bipartite pure states, it is known 
that both MLOCC and ELOCC are strictly more power- 
ful than ordinary LOCC [Hill. In the stochastic ver- 



sions of multiple-copy and entanglement-assisted trans- 
formations (SMLOCC and SELOCC, respectively) we 
are only concerned with non- vanishing success probabil- 
ity. For bipartite pure states, a transformation is realiz- 
able by SLOCC if and only if it is possible by SMLOCC or 
SELOCC, because of the multiplicativity of the Schmidt 
rank: Sch(* (g) $) = Sch(\E')Sch(<I>). Thus, there is no 
stochastic entanglement catalysis in bipartite systems. 

In this Letter, we advance both topics of multipar- 
tite tensor ranks and SMLOCC/SELOCC transforma- 
tions while demonstrating how results of the first have 
unexpected consequences for the second. As we show, 
since tensor rank is not multiplicative, there exist in- 
stances when the use of multiple copies or a catalyst can 
increase the conversion probability of some transforma- 
tion from zero to positive. In the first part of the paper, 
we describe the correspondence between homogeneous 
polynomials and symmetric states, and use it to bound 
the tensor rank of various multipartite symmetric states. 
In the second part, we derive some general properties of 
SMLOCC and SELOCC transformations and then use 
results from the first part to demonstrate the feasibility 
of certain SMLOCC and SELOOC transformations when 
their corresponding SLOCC conversions are impossible. 

Homogeneous polynomials and symmetric states. A sym- 
metric multipartite state is one that is invariant under 
any permutation of the parties, IW3) provides a tripar- 
tite example. For such a state we can ask not 
only about its tensor rank, but also about its symmetric 
tensor rank srk(';/'): the smallest number of symmetric 
product states {|'?!>i)^"}i=i,...,srk(?/;) to provide an expan- 
sion = E'^i'^^ |</>i)®"- To estimate srk(V') (and thus 
rk('0)) we introduce a correspondence between symmet- 
ric states and homogeneous polynomials. 

A homogeneous polynomial h of order N in d variables 
xi, . . . ,Xd is a linear combination of monomials x^- = 
^1 ■ ■ ' ('with a multi-index j — ji . . .jd), i-c. it has 



the form h — h{xi 



,Xd) 



$Z flj JJ^ , where 



the sum extends over all multi-indices with X^iLi 3i = ^■ 
Every homogeneous polynomial has a symmetric decom- 
position h YH=i\P^,i^^ + ... + Pn,iXd)^, with the 
minimum number pr(/i) of power terms. We refer to 
this number as the polynomial rank of h. The computa- 



tion and estimation of polynomial rank is a much-studied 
problem in algebraic geometry [l2, 17 1. 

Now, introducing a computational basis . . . , jd)} 
of the d-dimensional local systems Jia , a monomial x^ is 
associated with the Dicke state defined as 



I^(J)>:= 



N 

id 



1/2 



® \d) 



where Pgym is the projection onto the Bosonic (fully 
symmetric) subspace, Psym = m X^ttsSn ^'^i svmv 
extending over all permutation operators U.,^ of the N 
systems. General homogeneous polynomials (symmetric 
states) are associated by linear extension of the above 
since monomials (Dicke states) form a basis for the ho- 
mogeneous polynomials (symmetric states). That is: 



Observation 2 Every symmetric state 



d\®N 



uniquely associated with a homogeneous polynomial h(ip) 
of order N in d variables, and vice versa each homoge- 
neous polynomial h is associated with a symmetric state 
\h), such that h{D{j)) = x- and \x-) ~ \D{j)). Under 
this identification, symmetric tensor rank and polynomial 
rank are identical: pr(/i) = srk(/i). □ 

E.g., two copies of iWs) read | Wa)^^ ^ (|oo3) + 1030) + 
|300)) + (|012) + |021) + |102) + |120) + |201) + |210)), 
which is a sum of two Dicke states having correspond- 
ing homogenous polynomials 0:0X0X3 and xoXiX2- These 
have symmetric expansions xoXoXa = ^{{xo+xs)^~{xo — 
x^)^ - 2x1) E^nd X0X1X2 = ^((xo + xi+ X2)^ - {-xo + 
xi -\- X2Y — {xo — xi -\- X2Y — {xq -\- xi — X2Y), thus 
rV{Wf'^) < srk(VK3®2) < 7^ ^^^^^ jg ^jgi^^ ^ 

Using Observation [1] wc prove the following relations 
between unrestricted and symmetric tensor ranks. 

Theorem 3 (a) For multiqubit Dicke states 
\D{m,n)) := Ps„„(|0®™, 1®")) with m > n, 
rk(Z?(m, n)) = srk{D{m, n)) = m + 1, 

(b) for any N-partite symmetric state \4>), 
rk(V') < srk(^ ) < 2^-ir k(V>), 

(c) hm„^oo Aysrk(V'®") = hm„^oo 

Proof (a) The second equality follows from [13, 
Cor. 4.5] and it always holds that rk(Z)(m,n)) < 
srk(Z?(m, n)). So to prove the first equality, it suffices to 
show that the lower bound of rk(£'(m, n)) equals m + 1 
too. We use induction on n. For n = 1, the claim is 
true [l3l |. and we assume it holds for n — 1. Ignoring 
normalization, we can rewrite the state as \D{m,n)) = 
\D{m,n - 2))|11) + \D{m - l,n - 1))(|01) + |10)) + 
\D{m — 2,n))|00). Now we perform the global operation 
|1)(11| + i|0)((01| + (10|) on the last two systems which 
cannot increase the rank. The resulting (m + 71 — 1)- 
partite state is just the Dicke state \D{m, n — 1)) and so 
Tk{D{m, n)) > rk(D(m, n - 1)) = ?7i + 1. 

(b) Suppose that \Tp) has an optimal product state ex- 



pansion \A. 



\Ni). As lip) is symmetric, we 
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have 1^) = Psyn. {\A,) ® ... ® \N,)). But this is 

just a sum of rk(-0) Dicke states, each one correspond- 
ing to the monomial XAi....XNi. From [TtI . Prop. 11.6], 
pr(a:^;....a;Ar.) < 2^^^ which proves the claim. 

Part (c) follows directly from (b). □ 

Three copies of IW3). By Observation [21 the homoge- 
neous polynomial h{Wf^) can be written as |(a;oa;iX6 + 
X[)X2X5+xoX'iXi + xiX2X4)-\-^x'^XT . To computc its poly- 
nomial rank, we perform the following linear transforma- 
tions which do not change the polynomial rank: yi = 

Xi + X2 ~ X4, y2 ^ Xi - X2 + X4, m ^ -Xi + X2 + Xi, Z3 = 

l/2{x3 + X5), Z5 = l/2{x3 + xq), ze = l/2{x5 + xq). By 
using the fact that the polynomial rank is invariant under 
scalar multiplication, we can remove constant coefficients 
and obtain pT{h{W^^)) < pr^xoyize - yf) +w{xoy2Z^ - 
2/2) + pr(a;oy4Z3 - vl) + pr((yi + 2/2 + Vif + ajgxr) < 16. 
Here, the inequalities follow from [17|: Table 2]. With the 
lower bound Tk{Wf^) > 15 [11], we have 

Theorem 4 (a) rk(W3^^) = 15 or 16, 

(b) lim„_,oo \[rki^f^) < ^ « 2.52. □ 

This improves the previously best bound of T'k{Wf^) < 
21 [T^j. In particular. Theorem |4] implies that two tri- 
partite GHZ-type states with tensor rank 4 are sufficient 
to prepare three IVF3) states under SLOCC. 

Upper bound on the tensor rank of \Wn)'^"' . The A^- 
partite W state is defined as the Dicke state \Wn) = 

^(|o...oi)-H...-h|io...o)) e (c^)®^. AsIW^at)®" win 

be a linear combination of Dicke states, we can obtain an 
upper bound for rk(W^") by adding up the tensor ranks 
of each component Dicke state. Now each one of these 
corresponds exactly to a different way of separating n 
distinct excitations |1) into k = 1, 2, . . . , local states. 
This number is equal to the Stirling number of the second 
kind, namely S{n,k) = n ELo(-I)''"' ^^^re we 
only have to consider k < n as any larger number of 
parties is taken care of by the symmetrization. 

For example, S{3, 2) = 3 which implies that 
there are three ways of separating three excitations 
into two local systems, namely jO'^^, . . . , 0®^, 100, Oil), 
|0®3,...,0^3^010,101) and |0®3 . . . , 0^^^ 001, 110). By 
permuting the local states, each of these generates a 
Dicke state with corresponding monomials XQ~'^XiX3, 
Xq X2X5 and Xq XiXe, respectively. 

Since each of the S{n, k) monomials representing the 
same separation [n — > k) are related by a simple 
change in variables, each will have the same polyno- 
mial rank. Then by adding up all separations we ob- 
tain MWT) < E"="/'^'"^5(n,fc)pr(x^-^xi---Xfe) < 
^mm{w,n} ^^^^ fc) (i+niax{7V-fc, /t})2'=-i, where the sec- 
ond inequality follows from [13, Cor. 4.5 and Prop. 11.6]. 
In particular, this bound is of the form f(n)N + g(n) 
with some functions f{n) and g{n). In other words. 



Theorem 5 rk(W^") is upper bounded by a linear func- 
tion in N . Thus for large N , [Wn)^" can be prepared by 
LOCC from a GHZ-type state of rank linear in N . □ 

The large-n behavior of this bound is not very good, 
but based on a simple asymptotic consideration of the 
Stirling numbers for n « logA'^, we find that 

Corollary 6 hm„^^ s/MW^) < 0(log N). □ 

Lower bound on the tensor rank of \Wn)'^"' . 

Lemma 7 Any state of the form \n) = IWn-i)^'^ + 

ELi E.e5„ c.fcC/. (r^-i)«'=|o^_i)«("-^-)) 

SLOCC equivalent to |I^Ar_i)^". 

Proof We perform successively invertible SLOCC 
transformations on |f2), each transformation eliminat- 
ing a term in the double sum. For instance, apply- 
ing the transformation IM^at-i) — IM^at-i) — CTrk\ON-i), 
IOat-i) — >■ IOat-i) on |ri) by local invertible operators will 
eliminate the term C/^ (|iy7V-i)®''|0Ar_i)®(""'=)). The 
procedure is repeated on all terms in the sum until just 
IWat-i)'^" remains. □ 
To prove a lower bound, note that rk(W^") is the min- 
imum number of product states whose linear span con- 
tains the set S = {\Wn-i), IOat-i)}*"". Each of these 
product states can be substituted with an clement from 
S \ {\Wn-i)'^"} to yield a new set whose linear span 
also contains S. Thus, \Wn-i)'^"' is a linear combi- 
nation of elements from S \ {\Wn^i)'^"} and at most 
rk(iy^") - (2" - 1) product states. Thus by Lemma mwe 
get Tk{W^":^) < rk(VF®") - (2" - 1). As proven in Hi, 
for iV = 3, 2"+i - 1 < rk(W3®"). From these two inequal- 
ities, a simple inductive argument provides part (a) in the 
next theorem; part (b) then immediately follows after ob- 
serving that Theorem[5]reads rk(| Wjv)^^) < — 2 when 
n = 2. 

Theorem 8 (a) rk(W^®") > (iV - 1)2" - TV + 2, 

(h) Tk(W^^) = 3N -2. □ 

Multi-copy and catalytic SLOCC transformations. We 
now move on to the topic of SLOCC catalysis for 
multipartite entanglement transformations. Let H = 
^k=i ^fe ^^'^ = ®fc=i ^'k ri-partitc quantum sys- 
tems, and consider Hk and H'f. to be orthogonal to each 
other. Let jV'o) and be two vectors from T-L and H', 
respectively. Then the direct sum of lipo) and \(f)i) is 
given by |V'o)®|V'i) € H®H' C ^l^j^(Hk®H'f,). Notice 
that when = ®fe=iLfe|V'i) and \^2) = <»fc=iife|'/'2), 
we simply have |?!)i)®|02) = (8>Li(^fc®^A;)(l'/'i>® l'/'2)). 
By induction one can immediately show that the SLOCC 
ordering is preserved under direct sums. 

Lemma 9 If \tjjk) can be transformed into {(pk) via 
SLOCC, then \ipk) can also be transformed into 
0^^ \(f>k) via SLOCC. □ 
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We can use Lemma [9] to get a general relation be- 
tween SMLOCC and SELOCC. Assume that can 
be transformed into 10)*^" via SLOCC for some n > 1. 
Then by choosing |c) = 0)!^^ I^/))®"--'^® the resuh 

that (g) |c) can be transformed to (8) |c) via SLOCC 



follows from Lemma [5] So we get, similar to jl9l |2 



Theorem 10 // \^) can be transformed to \(f)) via 
SMLOCC, then the same transformation can also be 
achieved via SELOCC. □ 

By Observation [U to demonstrate the effect of entan- 
glement catalysis, wc only need to find a state \(f)) with 
the following property: rk((/)) = n and there is some fc > 1 
such that rk(0'^'^) < (n— l)*^. The source state l-ip) can be 
chosen as an n-partite GHZ state with tensor rank (n— 1). 
Such states do exist as proven in the previous section. 
In the following we shall provide two different construc- 
tions. The first class is given by the famous tripartite 
matrix multiplication tensor and the second one is given 
by the Wn states. By Theorem [10] these also suffice to 
show the existence of SELOCC transformations when the 
uncatalyzcd transformation is impossible. 

Theorem 11 Let ^ \^2)ab ® |$2)sc ® \<^>2)cA, 

where |$2) — |00) -f- and let \iP)abc be any general- 
ized CHZ-type state with tensor rank 6. Then the trans- 
formation of to l^''^') cannot be realized by SLOCC 
but can be realized by both SMLOCC and SELOCC. 

Proof It has been shown that 1$^'^)) is just the 2x2 
matrix multiplication tensor |21| . By a well known re- 
sult in algebraic complexity theory, rk( — 7 > 6 [2^ . 
Hence, cannot be SLOCC transformed into |(f>('^'). 
Now the best known algorithm for dx d matrix multipli- 
cation requires 0(d^ '^™) multiphcation steps . Hence 
the tensor rank of which corresponds to the al- 

gebraic complexity of 2" x 2" matrix multiplication, is 
0(22-376"). On the other hand, the tensor rank of l^-)"*" 
is simply 6" = 2^^°S2Q)'^ « 22-585"^ .^^hjch is larger than 
Q^2'^.376n^ for sufficiently large n. Thus we have con- 
firmed the existence of n (perhaps very large) such that 
rk(^'®") > rk(($(3))®")^ Y)oth SMLOCC and SE- 

LOCC arc possible. □ 
If we consider multipartite rather than tripartite state 
spaces, the W states provide much simpler examples. 

Theorem 12 For any N > 5 the transformation of 
\GHZ%^^) to \Wn) cannot be realized by SLOCC but 
can be achieved by both SMLOCC and SELOCC. Fur- 
thermore, two copies are sufficient in SMLOCC, and the 
catalyst in SELOCC can be chosen as \ Wki) ®\GHZ^~^) . 

Proof The result follows immediately from the facts 
that rk(iy^ 



N 



3iV - 2 in Theorem El and {N - if > 
2>N — 2 for iV > 5. One can easily see that the rank of the 
GHZ state can indeed be chosen as \\/'iN — 2] , which is 
much smaller than iV — 1 for 1 . □ 



Conclusions. We have shown that the theory of homoge- 
neous polynomials can be used to obtain insights on the 
symmetric tensor rank of symmetric states. Via this con- 
nection, we proved upper and lower bounds on the tensor 
rank for one and multiple copies of Wm states as well as 
the exact tensor and symmetric tensor rank of multiqubit 
Dicke states. We then proceeded to show that multi-copy 
and catalytic activation of otherwise impossible SLOCC 
transformations exists, using our results on Wn states to 
find explicit low-dimensional examples. 

Our work suggests several open questions which we 
leave for future investigation. First, given two states 
and 10), what are the necessary and sufiicient conditions 
such that |-0) can be converted to under SMLOCC 
and SELOCC? When |^) is a generalized GHZ state, 
the question becomes completely a matter of tensor rank 
multiplicativity. Asked in a different way, for some target 
state 10), when does there exist a state |-0) such that 
transformation jV') to |0) is possible under SMLOCC and 
SELOCC but impossible with just single copies. Another 
relevant problem is to determine the asymptotic tensor 
rank of IW3), and more generally of jW^v). Note that 
our lower bound of 2 coincides with the border rank 24 1 . 



It is conceivable that the asymptotic rank is 2 for all A^, 
but even an improvement of our logarithmic upper bound 
would be interesting. 
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